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Abstract. Chains of impact craters, or catenae, have been 
identified in Voyager images of Callisto and Ganymede. 
Although these resemble in some respects secondary crater 
chains, the source craters and basins for the catenae cannot be 
identified. The best explanation, proposed by Melosh and 
Schenk, is a phenomenon similar to that displayed by former 
comet P/Shoemaker-Levy 9: tidal (or other) breakup close to 
Jupiter followed by gradual orbital separation of the fragments 
and collision with a Galilean satellite on the outbound leg of 
the trajectory. Because the trajectories must pass close to 
Jupiter, this constrains the impact geometry (velocity and 
impact angle) of the individual fragments. For the dominant 
classes of impactors, short-period Jupiter-family comets and 
asteroids, velocities at Callisto and Ganymede are dominated by 
jovian gravity and a satellite's orbital motion, and are 
insensitive to the pre-fragmentation heliocentric velocity; 
velocities are insensitive to satellite gravity for all impactor 
classes. Complex crater shapes on Callisto and Ganymede are 
determined from Voyager images and Schmidt-Holsapple 
scaling is used to back out individual fragment masses. We 
find that comet fragment radii are generally less than -500 m 
(for ice densities), but can be larger. These estimates can be 
compared with those for the Shoemaker-Levy 9 impactors. 

Introduction 

Understanding the Shoemaker-Levy 9 impacts requires 
knowledge of individual fragment masses and densities, as these 
critical parameters governed the spatial and temporal deposition 
of energy and cometary materials in the jovian atmosphere. 
Pre-encounter diameter predictions ranged from up to several 
km for some fragments [Weaver et ai, 1994] to values an 
order of magnitude less [Scotti and Melosh , 1993]. Analyses 
of the events themselves have not yielded a consensus as of 
this writing. We argue here, however, that there is a historical 
record of similar comet fragmentations, expressed as otherwise 
enigmatic chains of impact craters on the Galilean satellites 
Callisto and Ganymede. Moreover, the trajectories that the 
crater-chain-forming fragments took are fairly well constrained, 
so by using modern crater scaling relations a reasonably 
definitive estimate of the masses of historical comet fragments 
can be obtained. We find that "large" fragments, defined here to 
be >1 km in diameter (assuming an ice density), are not ruled 
out by the historical record, but would be unusual. Below we 
discuss the geological record left by fragmented comets, the 
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crater shape and scaling issues that affect interpretation of this 
record, the probable fragment trajectories and velocities, and 
finally, our inversion for the fragment (and parent comet) 
masses and equivalent radii, and implications. 

Crater Chains 

Distinctive crater chains (or catenae) have been identified on 
Callisto, and subsequently, Ganymede. They are remarkably 
linear and unrelated to underlying geology (Figure 1 ; Schenk 
et ai , manuscript in preparation, 1995). At least 8 such 
chains have been identified on Callisto, 3 on Ganymede. 
Several linear chain-like features also exist, but at resolutions 
too poor to classify with confidence. Catenae craters on 
Callisto have well-developed raised rims and central peaks (or 
pits) (Fig. la), and two of the three recently detected chains on 
Ganymede [Melosh and Schenk , 1993] show unequivocal 
ejecta deposits, all strong evidence for an impact origin. 
Supporting evidence includes possible associated ejecta on 
Callisto (ejecta being more difficult to discern on Callisto’s 
less ice-rich surface) (Fig. la), chains that cross preexisting 
structures or structural boundaries ((i.e., basin rings, cratered 
terrain/grooved terrain) uninterrupted or influenced, and 
overlapping craters within several catenae (e.g., Figure lb). 

Among known impact features, the catenae best resemble 
secondary crater chains such as those radiating from large lunar 
basins [Wilhelms , 1987] or the Gilgamesh basin on 
Ganymede. Passey and Shoemaker [1982] originally judged 
the chains on Callisto to be secondary craters of collapsed or 
relaxed basins that had disappeared from the geological record, 
despite the high stratigraphic positions of some of the chains, 
but the search for possible source craters and basins did not 
meet with success [ McKinnon , 1983]. Fragments of 
impactors striking one satellite at grazing incidence and then 
traveling on to the next have also been suggested [Halfen et 
ai , 1990]. Whether or not this explanation works physically, 
the geometric cross-section ricochet trajectories must thread 


Figure 1. Voyager images of two crater chains on Callisto; 
same scale, north up. (a) 25-crater chain on the northern 
portion of the Valhalla mulliringed system, centered at 35° N, 
47° W (from FDS 16424.32, rectified); (b) 7-crater chain, 
centered at 49° N, 348° W (from FDS 16426.10, rectified). 
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through, and thus their probability, is substantially smaller 
than that appropriate to comets that break up near Jupiter and 
subsequently strike a satellite. 

The remarkable similarity of the Shoemaker- Levy 9 "string 
of pearls" led Melosh and Schenk [1993] to argue that the 
catenae on Callisto and Ganymede represent a fossil record of 
similar comet splitting events within the Jovian system, 
whether due to tidal stresses [e.g., Melosh and Schenk , 1993; 
Scotti and Melosh , 1993] or, possibly, to collisions with ring 
particles [Asphaug and Benz , 1994]. We adopt break-up near 
Jupiter as the initial condition for the formation of catenae. 

Crater Geometry and Scaling 

In order to determine the masses of the impactors that made 
the catenae craters, the mass excavated must be determined for 
each crater. Catenae craters range from 5 to 40 km in diameter 
(measured perpendicular to the catena axis when necessary to 
avoid overlap), and all (or nearly all) are morphologically 
complex, or collapsed. To determine the mass excavated, the 
transient, pre-collapse shape must be derived from the Final, 
observed shape for each. Shape cannot be measured reliably for 
catenae craters on Callisto or Ganymede. None are close 
enough to the terminator to use shadows, and strong albedo 
variations on Callisto make photoclinometric (shape from 
shading) slope estimates there unreliable [ Schenk , 1991]. 
Because of the strong similarity of catenae crater morphology 
to that of complex craters overall on Ganymede and Callisto, 
we use the general crater morphological statistics of Schenk 
[1991] to describe crater shapes. 

We model complex craters as "pie tins," with Hat floors 
(central peaks and pits have negligible volumes for our 
purposes), steep inner walls, and raised rims. Photoclinometric 
and shadow-length measurements of (normal) complex crater 
depths (d) and rim heights ( h ) on Ganymede from Schenk 
[1991] can be fit with power laws of the form d = 
0.26(±0.05)Df^ - 39 ± 0.06 anc j f x _ 0.09Df^ where D f is 
final crater diameter and all units are in km. The smaller set of 
shadow-length measurements of crater depths on Callisto are 
consistent with the Ganymede data, so because the geology and 
gravity on the two satellites are similar, we use the Ganymede 
fits for both. Inner rim slopes are taken at a constant 25° and 
outer rim flanks at 8°, an approximation, but one consistent 
with photoclinometric profiles [e.g., Schenk , 1993]. We as- 
sume volume conservation during the modification stage and 
neglect distributed rim subsidence, and model the transient ge- 
ometry as a paraboloid of revolution, extending the rim (at 8°) 
until a total depth/diameter of 1/2V^ [e.g.. Grieve and Garvin , 
1984] is achieved. This results in an approximate law relating 
transient to final crater diameters on Ganymede and Callisto: 

Df - 1.176D lr U08 (1) 

where D{ T is the transient diameter. This relationship is 
probably accurate to -15%, and the excavated volumes (and 
masses) to -50%. Equation (1) predicts that final craters in the 
5-40-km diameter range are -40-60% wider than their transient 
forms, substantially greater than the 30% enlargement assumed 
by Shoemaker and Wolfe [1982] in their treatise on cratering 
rates in the Jupiter system. Their enlargement factor was based 
on studies of collapsed lunar craters, but complex craters on 
Ganymede and Callisto are much flatter than lunar craters 
(which was not known in 1982), having collapsed to a greater 
degree [ Schenk , 1991, 1993]. The exponential dependence in 
equation (1) is, however, similar to that derived for lunar com- 
plex craters [see McKinnon et a/., 1991; Holsapple , 1993]. 

Once the transient crater volume V is determined, the 
impactor mass and radius a can be determined from scaling 


laws for a selected impact velocity u and impact angle from 
vertical 0, when u is replaced by ucosO [Chapman ana 
McKinnon , 1986]. Schmidt-Holsapple scaling is used (see 
Holsapple [1993]). Cratering efficiency (mass excavated 
divided by impactor mass) is determined as a function ol 
gravity-scaled size ;r 2 = gafu 2 for impacts into ice in both 
the strength and gravity regimes (g is surface gravity). The 
effective strength Y of ice is taken to be that used by 
Holsapple [1993] to model impacts into soft rock (see his 
figure 6), and is thus approximate; the surface density is 
assumed to be 1 g/crn^. Nearly all catenae craters form in the 
gravity regime, based on the range of fragment impactor sizes 
derived below and the impact velocities at Callisto and 
Ganymede (7T 2 lies between -3 x 10' 7 and 1 0' 5 ) ; hence the 
exact scaling in the strength regime is not so important. The 
absolute calibration of the cratering efficiency in the gravity 
regime remains uncertain, however, perhaps by a factor of 2 
(based on the range of previously published estimates). This 
dominates the overall uncertainty in our calculations. 

The Impactor Flux at Jupiter 

Impacts on Callisto and Ganymede are caused by comets and 
asteroids. Of these, short-period comets and asteroids on short- 
period-comet-like orbits (extinct comets) dominate the flux; 
Halley-family and long-period comets contribute perhaps 10% 
[Shoemaker et ai , 1986, Shoemaker , 1994]. The orbits ol 
97 short-period comets discovered through 1978, along with 
synthetic model orbits, were used by Shoemaker and Wolfe 
[1982] to determine a weighted encounter velocity « 4.5 
km/s at Jupiter's gravitational sphere of influence, only -0.35 
of Jupiter’s mean orbital speed u c j rc * 13 km/s. Gravitational 
focusing reduces this relative velocity for close encounters. 
Long-period comets, while rarer, encounter Jupiter at much 
higher speeds ~ V3t> c j rc ~ 22.5 km/s. New discoveries 
of extinct comets (asteroids) include higher inclination objects 
[ Shoemaker , 1994], which increase their weighted t;^. P R. 
Weissman [personal communication, 1994] recommends values 
closer to 0.6u c j rc = 8 km/s. At the risk of maximizing the 
derived impactor masses, we will adopt a value at the low end 
of the spectrum, = 0.3t) c j rc ~ 3.9 km/s, for short-period 
comets and dynamically similar asteroids or extinct comets. 
While the probability distribution of for short-period 
comets is uncertain, we shall see that it is not the dominant 
determinant of collision speeds on Callisto and Ganymede, and 
hence not the dominant source of uncertainty in our results. 

Fragment Trajectories 

The trajectories in Figure 2 illustrate encounters of short- 
period comets with = 0.3u c j rc that pass through Jupiter’s 
Roche zone, 1 to ~2.5(pj/<5) 3/3 /?j, where R j and pj are 
Jupiter's radius and density (= 1.33 g/cm 3 ), respectively, and <5 
is the cometary (or asteroidal) density, although tidal break-up 
may only occur within « 1 J(pjlS) 1/3 /? j [Sridhar and 
Tremaine , 1992]. Outbound the trajectories pass through 
Callisto's or Ganymede's position at either the ascending or 
descending node on the Jovian equatorial plane. Based on 2- 
body (Jupiter+comet) hyperbolic orbits, the angle the 
trajectories make with respect to the radius vector to Jupiter is 
constrained, -13° at Callisto and -18° at Ganymede for r min * 
1.5/fj; variations of several degrees are possible depending on 
r m in and the influence of Jupiter’s higher gravitational har- 
monics (/ 2 and J 4 ). With each satellite's motion (aberration) 
taken into account, projections of fragment impact trajectories 
shift somewhat onto the leading hemisphere of each satellite. 

Comet velocities (u c ) at each satellite's orbital distance are 
given by the root-sum-square of v « and the escape velocity 
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Figure 2. Trajectories in Jupiter’s equatorial plane for comet 
fragments with v ^ ~ 0.3u c j rc , from 3-body 

(Sun+Jupiter+comet) integrations of the equations of motion 
(courtesy L.A.M. Benner). Perijove distances are 1, 1.5 , and 
2.5 R } . The influence of J 2 and J 4 is illustrated and is 
maximized for equatorial plane trajectories. 


from Jupiter at the orbital distance (V2t> sat ), where u sat is 
the satellite's mean orbital speed; for the 2-body problem u c 
is independent of r m j n . Because fragment trajectories are 
nearly perpendicular to the satellite orbits, the approach 
velocity to a satellite is close to the root-sum-square of v c and 
u sa t. Further combined with a given satellite's escape 
velocity (u esc ), impact speed u of a fragment on 
Ganymede and Callisto is 

u = (3ur at + u* + u e 2 sc ) 1/2 (2) 

Escape velocity from Jupiter at Callisto's and Ganymede's 
radial distance is 1 1.6 and 15.4 km/s, respectively; t> esc < 3 
km/s for the satellites themselves. Hence, for v ^ = 0.3u c i rc 
~ 4 km/s, u is dominated by jovian gravity and satellite 
motion, and u = 15 and 19.5 km/s for Callisto and Ganymede, 
respectively. Even for = 0.6u c j rc , the depth of the Jovian 
gravitational well is such that u only changes by 10% at Cal- 
listo and 6% at Ganymede. Only for high-speed, long-period 
comets does v ^ become critically important, if not dominant. 

The wire frames in Figure 3 illustrate the locations of all 
the probable crater chains; five small chains in the original 
compilation of Melosh and Schenk [1993] are not included 
because their impact origin is equivocal. The projections are 
centered on the equator near 35° and 36° W, normal to the 
average of the projected trajectory vectors for Callisto and 
Ganymede, respectively (for r m \ n * 1.5/?j and v ~ 
0.3u c i rc ). All of the observed chains are located on the 
appropriate (tidally fixed) hemisphere of each satellite. The 
projected trajectory vector of an individual comet (its velocity 


relative to a satellite) depends on the inclination of that comet's 
trajectory with respect to the Jovian equator. As the 
inclination varies from prograde (direct) to retrograde, the 
position on a satellite’s surface where the trajectory vector is 
normal to the surface traces out an oval, for fixed r mm and 
V^. These ovals are shown in Figure 3 for r m j n = \.5Rj 
and = 0.3u c j rc . The direct and retrograde trajectories in 
Jupiter's equatorial plane are normal to a satellite's surface, and 
the cometary velocity minimized and maximized, at the 
easternmost and westernmost points on the ovals, respectively. 
Gravitational focusing of the trajectories by the satellites 
themselves is minimal (the cross-sections are altered by only 
1-2%) and is neglected in this calculation and below. 

The inclination of a trajectory also determines the 
orientation of a chain on the satellite, with two possible 
inclinations for each orientation. Impact sequence is the key to 
choosing among the two (see Figure lb), but this cannot as yet 
be done for most of the chains in Figure 3. Thus, for this 
letter we simply assume that u is given by the mean values 
above, and that 6 for each chain crater is the angular distance 
from the crater to the center of the appropriate oval in Figure 3. 
This introduces an uncertainty in u of -10% and in 6 of -10° 
to 15° for catenae on Callisto and Ganymede, respectively. 
Impact angle enters the scaling through cos0, however, so 
when 0 is small its uncertainty is of minor consequence. 

Fragment Masses and Radii 

The comet fragment masses as determined from the crater 
scaling and trajectory estimates above are shown in Figure 4 
with the Callisto and Ganymede data identified. The Callisto 
and Ganymede contributions are mutually consistent. 
Fragment mass estimates for all observed catenae craters range 
over 4 orders of magnitude, with a concentration near the center 
of the distribution. In contrast, mass estimates for individual 
catenae are limited in range to approx. I order of magnitude, 
and some are more nearly equal. Specifically, catenae whose 
number of craters best compares with the Shoemaker-Levy 
fragment train have a mass range of an order of magnitude. 
Thus, fragment masses vary from comet to comet, but individ- 
ual comets tend to split into more similar mass fragments. 
This would be consistent with tidal breakup models in which 
gravitational instabilities cause roughly equal mass accumula- 
tions [Asphaug and Benz , 1994; Solem , 1994], but not with 
theories that predict a characteristic size or size-spectrum of 
comet building blocks [Weidenschilling, 1994], unless the 
characteristic size(s) are much smaller than the fragments. 
While Gipul Catena and a shorter chain on Ganymede are the 
sole contributors to the high end of the mass distribution, 
without them the total range in fragment mass is still 10 3 . 



Figure 3. View of Callisto and Ganymede centered on the average trajectory vector for v ^ ~ 0.3u c j rc and 
r min = L5 R j. Crater chains are shown, as are loci of trajectory normal vectors (see text) and imaging limits. 
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Figure 4. Histograms of comet fragment masses and radii 
and (inset) parent comet radii. Radii assume a density of 1 
g/cm 3 . "Average" comet radius from Weissman [1991] is for 
comets brighter than absolute magnitude 1 1. 


Equivalent fragment radii in Figure 4 are calculated 
assuming an ice density (1 g/cm 3 ). Most fragment radii are 
small, -100-700 m, and less than many estimates of the 
Shoemaker- Levy 9 fragments. The distribution does have an 
extended tail to larger sizes, however, giving a median of -400 
m and a mean of -600 ± 500 m. The low end may be 
undersampled, but the resolution on substantial portions of the 
satellites is good enough to pick out smaller catenae, if they 
were abundant. 

Fragments larger than 1 km radius are found only for the 
largest catena on each satellite. If the Shoemaker-Levy 9 
fragments were large (>1 km diameter), the event would be 
historically unusual. Of course, noticing the comet may itself 
illustrate a selection effect; SL-9 may have been one of the 
historically larger, more massive comets. We note that 
adopting a different density for the fragments, e.g., 0.6 g/cm 3 , 
increases fragment size by a factor of 1.2, and if the satellite 
surfaces are denser than ice, fragment masses go up nearly in 
proportion. The total (mostly systematic) uncertainty in 
fragment mass may be a factor of -3, or -1.4 in radius. 

Reconstructed radii for the parent comets are also shown in 
Figure 4. Few comet radii and no comet masses have been 
measured directly, but our estimates overlap those in 
Weissman [1991]. Future work will involve more detailed 
examination of split comet trajectories and catena geology 
( Schenk et aL, manuscript in preparation, 1995), and 
beginning in 1996, new imagery from the Galileo Orbiter. 
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